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Abstract We study finite-time future singularities in J?(R,G) -gravity, where R and 
G are the Ricci scalar and the Gauss-Bonnet invariant, respectively. In particular, we 
reconstruct the F(G)-gravity and J^(R, G) -gravity models realizing the finite-time 
future singularities. We discuss a possible way to cure the finite-time future singular- 
ities in JP(R, G) -gravity by taking into account higher-order curvature corrections 
or effects of viscous fluids. 



1 Introduction 

Among the possible alternatives in order to explain the Dark Energy Issue, are the 
so called Modified Theories of Gravity Q] [2]. 

We would like to consider modified J^(R, G) -gravity, where the action is de- 
scribed by a function of the Ricci scalar R and the Gauss-Bonnet invariant G = 
R 2 -4R flv R^ v +R^ a R^ a . 

Many of modified gravity models bring the future universe evolution to finite- 
time singularities. Some of these singularities are softer than other and not all physi- 
cal quantities (scale factor, effective energy density and pressure) necessarly diverge 
at this finite future time. Note that singular solutions correspond to accelerated uni- 
verse, and often appear as the final evolution of unstable de Sitter space. 

The presence of finite-time singularities may cause serious problems in the black 
holes or stellar astrophysics J3). Thus, it is of some interest to explore the JP(R,G)- 
gravity models realizing singularities and if any natural scenario to cure such singu- 
larities exists. 

We use units of kg = c = h = 1 and denote the gravitational constant StiGn by 
K 2 = 87T/M P1 2 with the Planck mass M Pl = G N l/1 = 1.2 x 10 19 GeV. 
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2 The Model 



The action of (R, G) -gravity is given by 

J?(R,G) 



S = / <PxJ=t 



2k 2 



(i) 



where g is the determinant of the metric tensor g^ v and Jz? m atter is the matter La- 
grangian. The spatially-flat FRW space-time is described by the metric 



ds 2 = -dt 2 + a 2 (t)dx 2 



(2) 



where a{t) is the scale factor of the universe. 

From the action in Eq. (Q~|i, the FRW-equations of motion (EOM) are derived as 



3 2 

Peff = —nH , PeS : 



(2H + 3H 2 ) , 



(3) 



where p e ff and p e g are the effective energy density and pressure of the universe, 
respectively, and these are defined as 

Peff = --^ {p + ^2 i(^R R -&)- 6H ^R + G ^G - 2AH 3 # G ] J , (4) 



- [,%R AH&' R + 2&' R - G&' G 



16H (H + H 2 ) + 8// 2 Jv; 



(5) 



Here, H = d(t)/a(t) is the Hubble parameter and the dot denotes the time derivative. 
p and p are the energy density and pressure of matter, whereas &' R = dR^(R,G) and 
jFg = dQ^(R,G). For general relativity with &(R,G) = R, p e ff = p and p e ff = p 
and therefore Eqs. and (O are the Friedmann equations. Consequently, Eqs. (HJi 
and (0 imply that the contribution of modified gravity can formally be included in 
the effective energy density and pressure of the universe. 



3 Finite-time future singularities 

We consider the case in which the Hubble parameter is expressed as 

where h, to and Hq are positive constants and t < to because it should be for ex- 
panding universe. j3 is a positive constant or a negative non-integer number, so that, 
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when t is close to to, H or some derivative of H and therefore the curvature become 
singular. 

Such choice of Hubble parameter corresponds to accelerated universe, becouse if 
Eq.© is a solution of the EOM ©, it is easy to see that the strong energy condition 
(Peff + 3/? e ff > 0) is always violated when > 0, or is violated for small value of / 
when < 0. It means that in any case the singularity could emerge as final evolution 
of accelerated universe. 

The finite-time future singularities can be classified in the following wayJfQ: 

• Type I and Big Rip. It corresponds to /3 > 1 and j3 = 1. H and R (~ H 2 ) diverge. 

• Type II (suddenJS]). It corresponds to —1 < j3 < 0. R (~ H) diverges. 

• Type III. It corresponds to < j3 < 1. H and R (~ H) diverge. 

• Type IV. It corresponds to /3 < — 1 but j3 is not any integer number. Some deriva- 
tive of H and therefore the curvature becomes singular. 

We note that in the present paper, we call singularities for j3 = 1 and those for j3 > 1 
as the "Big Rip" singularities and the "Type I" singularities, respectively. 



4 Reconstruction method 

In order to study the finite-time singularities in J^(/?,G)-gravity, we use the recon- 
struction method [6]. 

We assume that the contribute of ordinary matter and radiation in expanding 
singular universe is too small with respect to the modified gravity, and we study 
the pure gravitational action of J£"(/?,G)-gravity, i.e., the action in Eq. © without 
^matter- In this case, it follows from Eqs. © and © that the EOM of &{R,G)- 
gravity are given byQ: 

24H 3 ,^' G + 6H 2 ^' R + 6H£' R + R&' R - G&' G ) = , (7) 
8// 2 #c + + 4H£' R + \6H£' G (H + H 2 ) 

+,^ R (4H + 6H 2 ) + ^-R^ R -G^ G ^0. (8) 

In the case of pure gravity, these two equations are linearly dependents. 
Moreover, we have 

R = 6 (2H 2 +H) , G = 2AH 2 (H 2 +H) . (9) 

It is easy to see that, in the case of Type I, II and III singularities, G and R tend to 
infinitive when t — > to in Eq.©, and in the case of Type IV singularities tend to zero. 

By using proper functions Z(t), P(t) and Q(t) of a scalar field which is identified 
with the cosmic time t , we can rewrite the action in Eq.© without ^f m atter to 



5= i J d 4 x^—g{Z{t)R + P{t)G + Q{t)) 



(10) 
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By the variation with respect to t , we obtain 

Z'(t)R+P'{t)G + Q'{t)=0, 



(11) 



from which in principle it is possible to find t = t(R,G). Here, the prime denotes 
differentiation with respect to t. By substituting t = t(R,G) into Eq. dTOb . we find 
the action in terms of J^(R,G) 



&(R,G) = Z(R,G)R + P(R,G)G + Q(R,G). 
We describe the scale factor as 

a(t) = a exp(g(t)) , 



(12) 



(13) 



where oq is a constant and g(t) is a function of time. By using the Eqs.(|7j>, $Q and 
( fT3l . the matter conservation law p + 3(p + p) and then neglecting the contribution 
from matter, we get the differential equation 

Z"(t ) + 4g 2 {t )P"(t) - g(t)Z'(t) + (8gg - 4g 3 (t))pi(t) + 2g(t)Z(t) = , (14) 

where the Hubble parameter is H{t) = g{t). By using Eq. (0, Q{t) becomes 



Q(t) = -2Ag\t)P'{t) - 6g 2 {t)Z(t) - 6g(t)Z'(t) . 



(15) 



It means that in principle, by solving Eq.(fl~4b on the singular solution (O, it is pos- 
sible to reconstruct ^(R, G) producing finite-time future singularities. 



In general, if Z(f ) ^ 0, &{R, G) can be written in the following form: 
&(R,G)=Rg(R,G)+f(R,G), 



(16) 



where g(R, G) ^ and f(R, G) are generic functions of R and G. From Eqs. (|7]i and 
([8]), we obtain 



Peff 



Peg 



2K*g(R,G) 



2AH i ^' G + 6H £ [R 



dg( R ' G ) , df(R,G) 



dR 



dR 



2K 2 g(R,G) 



+{&-R&' R -G&' G ) , 
%H 2 ^' G + 2&' R + AH&' R + \6H&' G (H + H 2 ) 



dR 



dR J 



)+&-R&' R -G&' c 



(17) 



(18) 



where p e ff and p e ff are given by the expressions in Eqs.©. 

The modification of gravity could be included into the Equation of State (EoS) 
of an inhomogeneus dark fluid with energy density p e ff and pressure p e ff 
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p eff = a>p eff + 9(H,A...), (19) 

where co is the constant EoS parameter of matter and &(H,H...) is a viscosity term 
given by 

Sf(ff,#...) = — 7 ) D -A {l + (o){&-R& R -G&' G ) 



2K 2 g(R,G) 

dg(R,G) df(R,G) 
dR dR 



[6H 2 (l + co)+4H] 



-H#g(4 + 6co) + m,<P' c [2H + H 2 (2 + 3(oj] + 2,£' R + 8// 2 #c j(20) 



The use of this equation requires that g(R, G) ^ on the solution. 
By combining the two equations in Eq.©, we obtain 

&(H,H...) = --L [2H + 3(1 +co)H 2 ] , (21) 



K 



where Eq.dT9b has been used. 



5 Singularities in J?(R, G) -gravity 
Big Rip singularity in F(G) -gravity 

As a simple example of reconstruction method, we examine the Big Rip singularity 
in Gauss-Bonnet F(G) -gravity, where ^(R,G) =R + F(G) and F(G) is a func- 
tion of Gauss-Bonnet invariant only. In this case, by putting Z(f) = 1, the action 
of Eq.dTOb can be written in terms of two proper functions P(t) and Q(t) and the 
variation with respect to t yields 

P'(t)G + Q'(t)=0, (22) 

from which we can find t = t(G) and the action in terms of R and F(G) 

F(G)=P(G)G + Q{G). (23) 

Eq.(fT4b and Eq.dTSb read 

4(^)-tfM^ + *)-o. 

Q(t) = -24 g '(t)^P--6g 2 (t). (25) 
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For the Big Rip singularity, j5 = 1 in Eq. ©. If we assume Ho = (the constant 
is negligible in the asyptotic singular limit t — > to), g(t) = hj (to — t) and the most 
general solution of Eq. d24l i is given by 



P[t) 



1 



-{2to-t)t + cx 



3-/i 



Ah (h -i; 

where c\ and c 2 are generic constants. From Eq. 



(to-t) 
3-h n 

, we get 



■C2, 



(26) 



Q(t) 



6h 2 



lAh 1 



('Q-t 
2h(h 



pi) -ci{t -t) 



2-h 



(to-ty 



(to- 



il!) 



Furthermore, from Eq. (f22b we obtain t in terms of G and, by solving Eq. d23l l, we 
find the most general form of F(G) which realizes the Big Rip singularity 



F(G) 



y/6/i 3 (l+/i) 

h(l-h) 



, — l+h 

VG + ciG^ +c 2 G. 



(28) 



This is an exact solution of the EOM in the case of Big Rip. The term cjG is a 
topological invariant. In general, if for large values of G, F(G) ~ OtG 1 / 2 , where 
a(^= 0) is a constant, the Big Rip singularity could appear for any value of h ^ 1. 
Note that C2G^ 1+/ '^ 4 is an invariant with respect to the Big Rip solution. 



Other types of singularities and more general &(R, G) -gravity case 

In a similar way, it is possible to reconstruct F(G) -gravity models in wich the other 
types of singularities could appear, when j3 7^ 1 in Eq.© and the scale factor, when 
Hq = 0, behaves as 

~h(t -ty- ' 



a{t) = exp 



0-1 



(29) 



We give some results. 

The asymptotic solution (in the limit t — > to) of F(G) when /3 > 1 is expressed as 



F(G) 



T2i 



(30) 



Hence, if for large values of G, F(G) ~ — a\/G with a > 0, a Type I singularity 
could appear. 

When j3 < 1, the asymptotic solution of F(G) becomes 



F(G)~a\G\?, y = 



_2P_ 

3/3 + 1 



(31) 



Finite-time singularities in modified ,9-{R, G)-gravity and singularity avoidance 7 

where a is a constant. If for large values of G, F(G) has this form with < y < 1 /2, 
we find < j3 < 1 and a Type III singularity could emerge. If for G — > — °°, F(G) 
has the form in Eq. d3~TT ) with — °° < y < 0, we find — 1/3 < j3 <0 and a Type II 
(sudden) singularity could appear. Moreover, if for G — s- 0~, F(G) has the form 
in Eq. OTb with 1 < y < oo, we obtain — 1</3<— 1/3 and a Type II singularity 
could occur. If for G — >• _ , F(G) has the form in Eq. (l3TT l with 2/3<7<l,we 
obtain — °° < j3 < — 1 and a Type IV singularity could appear. We also require that 
7 ^ In I (3n — 1 ), where n is a natural number. 

As a consequence, a large class of realistic models of F(G) -gravity, which re- 
produce the current acceleration and the early-time inflation, could generate future 
time-singularities, as for example 0: 

All this models contain power functions of G and for some choices of parameters 
could produce singularities. 

With reconstruction method it is possible to derive also more general &(R,G)- 
models producing finite-time singularities. For example, in the model J?(R,G) = 
R — aG/R, where a is a positive constant, could appear the Type I singularity, 
whereas in the model ,^(R) = R + aRY, where a and y are constants, could ap- 
pear Types II, III or IV singularities (for a review, see Ref.|6|). 



6 Curing the finite-time future singularities 

We discuss a possible way to cure the finite-time future singularities in F(G)-gravity 
and &(R, G)-gravity. In the case of large curvature, the quantum effects become im- 
portant and lead to higher-order curvature corrections |9|. It is therefore interesting 
to resolve the finite-time future singularities with some power function of G or R. 

We consider the description of modified gravity as inhomogeneous fluid. If some 
singularities occur, Eq. d2TT i behaves as 

TO «...),-ffi^(, -,)-» + ^h-,r'-'. <33» 

One way to prevent a singularity appearing could be that the function @(H,H...) 
becomes inconsistent with the behavior of Eq. ( [331 in the singular limit (t — s- to). 

Let us consider a simple example in order to cure Big Rip singularity in F(G)- 
gravity. Suppose that for large values of G, 

R + F(G^°°) — >R + yG m , m^l, (34) 

with 7^0. For H = h/(to — t), namely the Big Rip case, we have 
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a 



(35) 



(to - 1) 



Hence, if m > 1/2, &(H,H...) tends to infinity faster than Eq. (|33l and we avoid 
this kind of singularity. 

As general results, we find that the term yG" with m > 1/2 and m ^ 1 cure the 
singularities occurring when G ±°° (Type I, II and III). Moreover, the term yG m 
with m < cure the singularities occurring when G^O" (Type II, IV). 

In f(R) -gravity (namely, R plus a function of R), by using the term yR m , the 
same consequences are found. The term yR m with m > 1 cures the Type I, II and III 
singularities. The term yR m with m < 2 cures the Type IV singularity. 

Within the framework of &(R, G)-gravity, we can use the terms such as G m /R" 
to cure the singularities. For example, we can avoid the Type I singularities if the 
asymptotic behavior of the model is given by yG'"R", with m,n > 0. 



7 Effects of viscous fluid in singular universe 

As the last point, we explore the role of perfect/viscous fluids within singular modi- 
fied gravity, investigating how the singularities may change or disappear, due to the 
contribution of quintessence or phantom fluids. 

We consider the class of modified gravity J^(R,G) = R + f(R,G), where f(R, G) 
is a function of the Ricci scalar R and the Gauss-Bonnet invariant G, and we suppose 
the presence in the universe of cosmic viscous fluid, whose EoS is given by 



where p and p are the pressure and energy density of fluid, respectively, and ft) is the 
EoS parameter. £ is the bulk viscosity and in general it could depend on p, but we 
will consider the simplest case of constant viscosity only (for more general cases, 
see Ref. lfTOl . ifTTI ). On thermodynamical grounds, in order to have the positive sign 
of the entropy change in an irreversible process, £ has to be a positive quantity. 
The FRW-equations of motion are: 



The modified gravity is formally included into the modified energy density Pq and 
the modified pressure pq, which correspond to Eqs.dT7T>-dT8b for g(R, G) = 1 . 
The fluid energy conservation law is a consequence of the EOM ( l37l i: 



p = (op-3HC, 



(36) 



Pg + P 




(37) 



p+3//p(l + ft)) = 9// 2 C. 



(38) 



The presence of fluid could influence the behaviour of singular /(7?,G)-models 
(i.e. models that in absence of fluids produce some singularities). We will check the 
solutions of the fluid energy density when H is singular. 
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In the non-viscous case £ = (perfect fluid), the solution of Eq.d38t assumes the 
classical form: 

p=p«a{t)-^ +m \ (39) 

where po is a positive constant and a(t) is the scale factor of the universe. By com- 
bining Eq.([39l with Eq.(|29t, it is easy to see that for j3 > 1 (Type I singularity), p 
grows up and diverges exponentially if (0 < — 1 . In the presence of phantom fluid, 
the EOM d37l ) become inconsistent with respect to the singular form of Hubble pa- 
rameter in Eq.©, and the Type I singularity is not realized in f(R, G)-gravity. 

When < /J < 1, the fluid energy density p is avoidable on the Type III singu- 
lar solution, whereas for Type II and IV singular models (j3 < 0), the presence of 
quintessence or phantom fluids makes the singularities worse. In particular, in the 
case of P < — 1, the dynamical behaviour of Eqs.(l37T> could become inconsistent, 
because p behaves as (fo — f ) and it is larger than the time-dependent part of H 2 
(~ (to-t)~P) when f fo- 



Constant viscosity 

Suppose to have the bulk viscosity equal to a constant, £ = £o, and the Hubble 
parameter in the general form of Eq.©. The asymptotic solutions of Eq.(f38b in the 
singular limit f — s- fo are: 



3/1 Co 



r (l + co)(t Q -t)P ' 

%oh 2 
(2j3-l)(f -f) 2 ' 3 - 1 ' 

%//oCo , 3i/ Co 



(p-l)(t -t)P- 1 1 + 63 



P>1, (40) 

1>J3>0, (41) 
0>p\// ^0. (42) 



In the first and second cases (j3 > 0), it is possible to see that p diverges more slowly 
than H 2 , so that viscous fluid does not influence the asymptotically behaviour of 
Types I and III singular models in Eqs.d37li. due to the constant viscosity. 

In the third case, we consider fluids that tend to a non-negligible energy density 
when J3 < 0. It automatically leads to Hq ^ in Eq.© and p behaves as in Eq.(l42l). 
Large bulk viscosity Co becomes relevant in the EOM. Moreover, if CO < — 1, the 
effective energy density (namely, Pg + P) could be negative and avoid the Type II 
and IV singularities for expanding universe (where Hq > 0). 
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8 Conclusion 

We have investigated the finite-time future singularities in F(G)-gravity and J?(R,G)- 
gravity. We can reconstruct the F(G)-gravity and JP(R, G) -gravity models in which 
the singularities could appear. Note that all types of future-time singularities could 
appear in J^(/?,G)-modified gravity. In addition, we have discussed a possible way 
to resolve the finite-time future singularities in F(G) -gravity and J?(R,G) -gravity 
under quantum effects of higher-order curvature corrections or the presence of per- 
fect/viscous fluid in the universe. 

Acknowledgements I wish to thank K. Bamba, S. D. Odintsov and S. Zerbini, whom the paper 
'Finite-time future singularities in modified Gauss-Bonnet and ,i?(R,G) gravity and singularity 
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